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ABSTRACT 


SMITH, DEWEY M. A Comparison' of Experimental Heat Transfer Coef- 
ficients in a Nozzle With Analytical Predictions From Bartz'-s 
Methods for Various Combustion Chamber Pressures in -a Solid 
Propellant Rocket Motor. (Under the direction of JOHN NOBLE 
PERKINS . ) 

The experimental heat transfer coefficients measured in the 
nozzle of a small solid propellant motor are compared to the predic- 
tions from D. R. Bartz's Nusselt number correlation equation and his 
technique of solving the boundary layer momentum and energy equations 
simultaneously for the heat transfer coefficient. The propellant was 
a composite of ammonium perchlorate- and polybutadiene acrylic acid 
and the average motor chamber pressures were 220, 410 , and psia. 
The nozzle was made with a steel casing and a. ZTA graphite throat 
insert. Measurement locations were at a local, to throat area ratio 
of I .785 in the convergent section, at the throat and’ at a local to 
throat area ratio of 2.569 in the divergent- section for the chamber 
pressures of 220 and 4l0 psia. The measurements were made at a 
local to throat area ratio of 2.067 in the convergent section, at the 
throat, and at a local to throat area ratio of 3 - 76 ^ in the divergent 
section for the chamber pressure of Jk2 psia. Test measurements 
consisted of temperature responses of five thermocouples mounted on a 

line perpendicular to the> heated' -surface of the nozzle at each of 

, » ♦ 1 

the measurement locations . The heating rates and. the corresponding 
convective heat transfer coefficients were determined by using the 
thermocouple data as input to- a finite difference heat balance -program. 



The results of the comparisons showed the experimental data from the 
convergent region and throat to he consistently lower than the predic- 
tions made using Bartz's two techniques. It was found that the experi- 
mental data in the divergent sections could ,be correlated by 

J ' . 

„ t 

evaluating the skin friction coefficient at the ’free stream temperature 
and using it in the simultaneous solution. 
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INTRODUCTION 


During recent years, attempts have been made to analytically 
predict heat transfer rates in the combustion chamber and along the 
nozzle wall of solid propellant rocket motors . The prediction tech- 
niques were developed by researchers in the liquid propellant rocket 
engine field and they have been fairly successful in defining the 
heating loads for this type engine. 

The two most widely used methods for predicting the heat transfer 
rates were set forth by D. R. Bartz. The first method Is a Nusselt 
number correlation equation which was developed before the advent of 
high speed computers and is still used today for rapid estimations. 

The second method solves the boundary layer momentum and energy equa- 
tions for the heat transfer coefficient and requires the use of a 
computer . 

These techniques are now being applied to solid propellant rocket 
motor nozzles. The experimental data on heat transfer rates for this 
type motor are not extensive and therefore it is difficult to determine 
the applicability of Bartz' s techniques to solid propellant motors. 

Most of the nozzles on solid propellant motors are not externally 
cooled and their design is based on the materials in the nozzle wall 
being able to absorb the heat transferred' from the exhaust gases. 

During the times Immediately after ignition, severe temperature 
gradients are set up through the nozzle wall and this condition must 
be taken into consideration in the design. Accurate prediction of the 
heating loads is desirable so that the nozzle may be designed 


efficiently. 
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To provide additional experimental data, the heating rates in a 
converging - diverging nozzle on a small solid propellant rocket motor 
■were measured at a point in the convergent region, at the throat, and 
at a point in the divergent region. The average chamber pressures were 
220, 4l0 , and 7^2 psia which represent the range of chamber pressures 
used in full scale motors. The heating rates were determined from 
the temperature responses of five thermocouples mounted on a line 
perpendicular to. the heated surface of the nozzle at each of the 
measurement locations. The experimental heat transfer coefficients 
were then compared with predictions from the two techniques of Bartz. 
This thesis presents these comparisons and the intention is to provide 
some basis for applying the Bartz techniques to solid propellant rocket 
motor nozzles. 
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REVIEW OE THE LITERATURE 

Bartz ( 4 ) has documented In detail the developments of both the 
analytical and the experimental research in rocket nozzle heat transfer. 
Published in 1965 , it provides a convenient catalog of .the most recent 
experimental results and techniques and was referred to often in the 
performance of the work for this thesis. 

The initial prediction method of rocket nozzle heat transfer was 
based on the turbulent pipe flow- heat transfer correlation equations 
of McAdams ( 13 ) ♦ This method assumed the flow in the nozzle to be 
fully developed with each point on the nozzle contour assumed 1 to be 
preceded by a long. pipe. Sibulkin and Bartz (15, 2 ) were the first, 
to treat the nozzle heat transfer problem with a boundary layer 
approach by making use of the integral momentum and energy equations 
(Sibulkin's was an incompressible flow analysis) . These treatments 
were an improvement over the McAdams equations since the flow in 

4 V * t 

’ * > 

rocket nozzles is. not usually fully developed. The main difference 

in Bartz 's initial boundary layer analysis and the one used today is 

the method of solving the boundary layer equations. The original 

analysis was done before the advent of high speed computers. Bartz, 

pointed out the equations were interdependent upon the ratios of wall 

T w 

temperature to stagnation temperature, — , and the temperature to 

■*•0 

velocity boundary layer thickness, ^ . However, to simplify the 
mechanics of solution of the equations, he assumed initial values for 
these ratios and solved the two equations separately. The boundary 
layer equations were reduced to linear ordinary differential equations 
with variable coefficients. 
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The other basic assumptions in the analysis were: 

(!) 1/7 power law profiles of both velocity and the difference 

between stagnation temperature and wall temperature in the turbulent 
boundary layer, 

(2) the local skin-friction coefficients along the nozzle are 
the same as those on a flat plate for the same boundary- layer thick- 
ness, and 

(3) Reynolds analogy between momentum transfer and heat transfer 

applies for the nozzle boundary-layer flow. Elliot et al. (8) 

developed a computer program in 1963 that solves the boundary-layer 

equations simultaneously by an iterative method which allows the 
T w A 

ratios, — and — , to vary along the nozzle wall. 

-i-o o 

A Eusselt number correlation equation was also developed by 
Bartz ( 3 ) for the purpose of making calculations of the local heat- 
transfer coefficients by hand. This method was based on the solution 
for the heat-transfer coefficient from the original boundary-layer 
•analysis which showed the local coefficient to be strongly dependent 
upon the local mass flow rate. The boundary- layer solutions also 
showed that the local diameter should be used as the characteristic 
length in the correlation equation. The proportionality constant was 
obtained by matching the heat-transfer coefficients at the throat of 
a particular nozzle with the coefficient that was calculated from the 
boundary- layer analysis. This correlation equation is still used 
today for rapid calculations of the local heat-transfer coefficient. 
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Numerous experiments have been performed to determine the accuracy 
of the methods. These experiments, for the most part, have involved 
heat transfer measurements in nozzles with heated air or various liquid 
propellant exhaust gases as the working fluid. The heated air experi- 
ments are unique in that the inlet conditions to the nozzle can be 
controlled thus eliminating such combustion effects as secondary flows, 
oscillations in pressure, or free stream turbulence that can occur in 
rocket motors. Of particular interest, are the data of Back, _et al. (l) 
which are from heated air experiments with stagnation pressures ranging 
from 30 to 250 psia and stagnation temperatures over the 1000 -2000° R 
range. The nozzle was a 30-degree half-angle convergent, 15-degree 
half-angle divergent nozzle which is similar to the nozzles used in the 
set of tests reported on here. 

The boundary layer analysis predicted the heat transfer coefficient 
accurately throughout the nozzle whereas the correlation equation was 
approximately 50 percent high for chamber pressures in the range 75 - 
250 psia. One other interesting point was the indication that the 
boundary layer turbulence seemed to decay back toward transition of 
laminar flow near the throat in the lower stagnation pressure tests 
and was reflected in the heat transfer results. 

The data obtained by Kolozsi (10), using air at stagnation pres- 
sures of 225 and 370 psia and stagnation temperatures of 1100 - 1200° R 
in a convergent-divergent nozzle, indicated the correlation equation 
predicted coefficients too high (^5 percent high at the throat) and 
the boundary layer equations were accurate throughout the nozzle. 

Fortini and Ehlers ( 9 ) found that both methods predicted the heat 
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transfer coefficients accurately in a Rao-design divergent section 
nozzle using air at a stagnation pressure of 300 psia and a stagnation 
temperature of l600° R. Their results also showed that two-dimensional 
flow must be considered in this type nozzle instead of one- dimensional 
flow that can be used in most convergent- divergent nozzles. 

Welsh and Witte (19) used a 1^0^ - hydrazine liquid propellant 
rocket engine to gather heat transfer data and compared it with only 
the correlation equation for stagnation pressures between 80 - 290 psia. 
They found the predictions to be considerably lower in the convergent 
region when compared with the experimentally determined coefficient, 
from 80 percent above to ^5 percent below in the throat region, with 
the best correlation in the divergent section. They theorized that 
the effects of combustion in the vicinity of the nozzle inlet 
influenced the flow in the convergent and throat region. Convergent- 
divergent nozzle configurations were used in these tests with varying 
contraction ratios. 

Witte and Harper (20) used the same liquid propellant engine used 
by Welsh and Witte (19) with nozzles over an extended range of contrac- 
tion and supersonic area ratios. They had the same general results 
as in ( 19 ) • They also concluded that the wide variation of data in 
the throat region of ( 19 ) was due to the transitional tendencies of 
the turbulent boundary layer here. This tendency was felt to be 
caused by the acceleration of the flow. 

Lee (12) obtained experimental heat transfer data from a solid 
propellant motor with an uncooled molybdenum nozzle that was found to 
agree with Bartz's correlation equation. The data could be correlated 
also by assuming the skin friction coefficient to be dependent on the 


momentum thickness and numerically integrating the boundary layer 
momentum equation. Brinsmade and Desmon (5) conducted tests with a 
solid propellant motor at stagnation pressures between l6o - 300 psia 
and a stagnation temperature of ^900° R and found that the data at the 
throat could be correlated by using laminar heat transfer equations. 
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THEORY 

The integral momentum and energy equations of the turbulent 
boundary layer can be derived by writing momentum and energy balances 
on a control volume in the vicinity of the wall where there are viscous 
effects in a real fluid. This was the method used by Bartz (^4-) . The 
derivation is based on the definitions of displacement, momentum, and 
energy thicknesses as the deficiencies in mass, momentum, and energy 
caused by friction and heat transfer. These deficiencies are deter- 
mined by comparing real flow with potential flow near the wall where 
the mass flow rates are made equal in the control volume for these two 
types of flow (see Fig. l) . 

The basic definitions are: 

(l) Displacement thickness 



which is the difference in thickness in the two control volumes in 
order to have the mass flow rates equal. 

(2) Momentum thickness 

* \ 

which is the thickness of potential flow which has a momentum flux 
that is equal to the difference between the potential and real flow 
momentum fluxes for the same mass flux. 
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(5) Energy thickness 



(3) 


which is the thickness of potential flow that has an enthalpy flux 
equal to the difference between the enthalpy fluxes of the potential 
and real flows for the same mass flux . In forming these definitions, 
the time-mean flow density, pu, has been represented by the product of 
the mean values, p and u, and ignoring the cross- correlation terms. 
This can be done on the assumption the correlation terms cancel out 
when integrated over the boundary layer as suggested by Shapiro (l4) . 

Using these definitions in the momentum and energy balances on 
the control volumes in the potential and real flows, the following 


equations are formed: 

(l) Integral Momentum Equation 

5* 


d0 

dx 


0 , 


- 8 


1 + 


~S~ dU 1_ d(pu) 1 dr 

U dx pU dx r dx 


where is defined as 


2t 


C f “ 


w 


pU 2 


W 

(5) 


(2-) Integral Energy Equation 




H 

1 d(pu) 1 dr 1 v 
pTT ”cEx + r ctx “ T Q 35T 


(6) 


where C, is defined as 

h 


= * UC p( T av - 


(7) 
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The basic assumptions used by Bartz in solving these equations 
are as follows; 

(1) The flow is axi symmetric and steady, the forces acting on 
the gas are the pressure gradient and shin friction at the wall, and 
the boundary layer is small compared to the distance from the axis of 
symmetry. 

(2) The flow through the nozzle is reversible and adiabatic with 
the change in total enthalpy of the’ gas due to the heat flux’ to the 
wall . 

( 5 ) The gas is perfect, has a. constant Prandtl number, and its 
viscosity is related to the gas temperature raised to a. power, 

(k-) The shin- friction coefficient and the Stanton number are 
the same as they would be on a flat plate at the same free- stream 
conditions, wall temperature, and momentum thichness, 

(5) The Stanton number for unequal momentum and energy thiehnesses 
is that for equal thicknesses multiplied by ||*J , where (n) is an 

"interaction component." The Stanton number for equal momentum and 
energy thicknesses is related to the skin- friction coefficient by 
von Karman's form of Reynold's analogy 



(6) Heat transfer has either no effect on the skin- friction 
coefficient, and C-f is the same as for adiabatic flow, or has an 
effect and the is the same for adiabatic incompressible flow with 
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the density and viscosity evaluated at the arithmetic mean of the wall 
and free stream static temperatures. 

(7) The boundary layer velocity and temperature distributions 
are l/7-power profiles. 

Values for the skin friction coefficient are taken from data of 
adiabatic flow over flat plates correlated by Coles (7) * Coles found 
the data could be represented by one curve of Cf versus CpRg where 
the low speed value, is related to the actual C^ a by 



where T g is a temperature within the boundary layer which is found 



and m is the exponent in the viscosity relationship, ,p ~ T 131 . The 
R^- is related to R Q by 


c fa R e 



R S 



( 11 ) 


A film temperature correction _>can be made by evaluating the gas 
properties p and u, at the arithmetic mean temperature of the free 

► 4 

stream temperature and the wall temperature. When this is done the 
relationship between C^ and is 



( 12 ) 
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The term, — , in the momentum equation can be evaluated from the 
integrals -which define 8* and 0, equations (l) and (2), by using the 
assumed l/7-power profiles for the velocity and temperature ratio in 
the boundary layer. The limits on the integrals for these thicknesses 
as well as for the momentum thickness, 0, are expressed in terms of 8 
and A, which are not necessarily equal. Therefore the boundary layer 
momentum and energy equations are solved iteratively for 8 and 0 
in order to determine the proper ^ as well as the proper skin 
friction coefficient. 

The mass rate of flow per unit area through the nozzle can be 
determined from one dimensional isentropic relationships. By 
specifying the nozzle contour, wall temperature, and chamber condi- 
tions, the heat transfer coefficient at any point along the nozzle 
wall can be found from solving the boundary layer equations for the 
Stanton number, C : 

' Vi 


c h ' puc 


% 


pUfT 


7r - t T 

^ aw w J 


(15) 


In the solutions for h developed by Bartz (2), it was found 

8 

that the heat transfer coefficient was ,a strong function of the mass 
flow rate per unit area, pU-, From this, Bartz. developed a nondimen- 
sional equation in the form 


N. 


TTu 



(14) 


for the determination of h . The analysis by Bartz (2) showed that 

8 

(a = 0.8) and the exponent, b, was evaluated to be 0.4 from 
von Karman's modification of Reynold's analogy for Np r = 1. The 
characteristic length in the Reynolds number is the local diameter 
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which varies approximately with the boundary layer. One notable 
exception to this variation is in the entrance region of a nozzle 
where the boundary layer thickness may be small. 

The constant C was evaluated by determining the heat transfer 
coefficient at the throat of the nozzle shown in Figure 2, with 
the boundary layer analysis and solving for C in equation (l4) using 
this h . To insure that this equation would apply to other nozzle 
contours and conditions, a factor — , found from nozzle similarity 

VoJ 

studies in (2) was multiplied into the equation. If it is assumed 
that the specific heat, Cp, and the Prandtl number, Fr, are constant 
with temperature and that the gas properties p and q are evaluated 
at a reference temperature to account for compressibility and/or heat 
transfer effects, equation (l4) can be expressed as 



This equation can be used to obtain a rapid estimation of the 
heat transfer distribution. It has been found to give good results 
except in the entrance region of nozzles with thin boundary layer 


thicknesses . 
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DESCRIPTION OP EXEER3IIE1-IT.AL TECHNIQUE 

Discus sloa of Experimental Apparatus 

The rocket engine used -to obtain the experimental heat transfer 
data is shown in the schematic in Figure 3- The solid propellant was 
a composite of 83*3 percent by weight ammonium perchlorate (oxidizer), 
lk .k percent polybuta'diene acrylic acid (fuel) and 2. 3 percent 
stabilizer. This propellant was selected because there are only small 
amounts of solid 1 particles (carbon residue) in the exhaust products, 
which minimizes heat transfer by radiation, and because its combustion 
chamber temperature is relatively low (about k500° F) . There is very 
little dissociation of the exhaust gases in this temperature range. 

Average chamber pressures of 220, klO, and 7^2 psia were obtained 
by varying either the throat area or the volume of the combustion 
chamber and the amount of propellant. Table 1 .gives the pertinent 
parameters for 'each of the chamber pressures used in the experiments 1 . 

It was .desirable to have a constant chamber pressure over the data 
taking period. To approximate this condition, the propellant was 
cast in a hollow cylinder configuration with the inner cylindrical 
surface being used as the exposed, burning surface. Figures k, 5? and ■ 

6 show the pressure traces for- the three^ test conditions. It can be 
seen that a good approximation of' constant chamber pressure with time 
was achieved. 

The convergent -divergent nozzle, was made 'with' a ZTA graphite 

» 1 1 ' 

* 

insert and a steel housing. The. internal surface contour was made smooth 
and continuous which allowed ’an assumption to be made that the 
beginning pf the boundary layer coincided with' the beginning of 
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convergent section of the nozzle. The ZTA graphite was chosen as the 
material for the insert because of i.ts machinability, resistance to 
erosion, and. its thermal properties. Also by using the graphite 
insert, the experimental nozzle configuration was similar to nozzles 
used in full scale motors. 

Heat transfer measurements were made in the convergent region, at 
the throat, and in the divergent region of the nozzle for the three 
test conditions. A degree of redundancy was obtained by using two 
calorimeters diametrically opposed in the divergent region. These 
calorimeters also were used to determine whether the gas flow through 
the nozzle was concentric. The exact locations of all four calori- 
meters in the nozzle for each of the. tests, are shown in Figure 7* 

The calorimeters were made with a ZTA -graphite core and a, silica 
phenolic insulating sleeye. As can be seen in Figure 3 j the graphite 
core was exposed directly to the flow of gases. By making both the 
insert and the calorimeter core out of the same material, there was 
very little disruption of the temperature distribution along the nozzle 
wall. The insulating sleeve was used to direct the heat flow along the 
longitudinal axis of the calorimeter. This allows the assumption of 
one- dimensional heat transfer to be made in the data reduction. The 
graphite core of the calorimeters was instrumented with five thermo- 
couples. The thermocouple nearest the exposed surface was composed 
of tungsten - 5 percent rhenium and tungsten - 2 6 percent rhenium wires 
5 mils in diameter. The other four thermocouples were made of 
platinum and platinum -13 percent rhodium wires also 5 mils in diameter. 
Figure 8 shows the locations of the thermocouples along the 
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longitudinal axis of the core. All calorimeters used in the tests 
were similarly constructed. 

By using five thermocouples in each calorimeter, a temperature 
distribution through the calorimeter could be obtained even if there 
were random failures of individual thermocouples. Also the thermal 
diffusivity of the graphite could be checked by specifying the tempera- 
ture-time history of one thermocouple in the data reduction analysis 
and comparing the calculated temperature distribution through the 
calorimeter with the distribution obtained from the experiment. This 
procedure will be demonstrated in a later section. 

Ideally, the surface temperature of the calorimeter should be 
directly measured. However, due to the extreme thermal environment 
at the surface, this measurement is very difficult to make. The 
surface temperatures of the calorimeters were determined by locating 
a thermo couple as closely as possible to the surface and extrapolating 
the temperature data to the calorimeter exposed . surface. 

All of the heat transfer measurements were taken from static 
firings of the rocket engine. Figures 9 and ’10 show the engine in 
position on the thrust stand. The thermocouple data as well as chamber 
pressure and thrust data were recorded on magnetic tape with a 
computerized data acquisition system. 

Thermal and Physical Properties 

The exhaust gas constituents were determined- by assuming the 
products of combustion to be in chemical equilibrium at the prescribed 
pressure and enthalpy in the combustion chamber. Table 2 presents the 
conditions in the combustion chamber for the three tests and the 
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resulting mole fractions of each exhaust gas constituent as well as 
the specific heat ratio, j, and the specific heat at constant pres- 
sure, Cp, for the gas mixture. These values were determined with the 
use of a computer program described by Zeleznik and Gordon (21) . The 
gas mixture was assumed not tc> vary in the nozzle (frozen flow) and 
local conditions at the test measurement locations were determined 
by expanding the mixture isentropically through the nozzle. 

The viscosity and thermal conductivity of the gas mixture were 

determined by using the equations suggested by Brokaw-(6): 

V' 


fI MiJ 


=1 

1=1 
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h. 

ij \ 

d=l 
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( 15 ) 



(16) 


Sax “ Sax + Skx 


( 17 ) 


These equations express the mixture viscosity and thermal conduc- 


tivity in terms of the viscosity and thermal conductivity of the 
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mixture costituents. The constituents that are starred in Table 2 
were accounted for in these calculations since they represent the 
major quantities in the mixture. The viscosity and thermal conductivity 
of each component were taken from data presented by Svehla (l6) . Their 
values were evaluated at the combustion chamber temperature. 

The properties of the ZTA graphite used in the calorimeter core 
that are required are the density, constant pressure, specific heat, 
and the thermal conductivity. The density and specific heat data 
were taken from (17) • The thermal conductivity data were taken from 
the results of Wagner and Dauelsberg (l8) . Figures 11 and 12 show 
the variation of the specific heat and thermal conductivity of ZTA 
graphite with temperature. The thermal conductivity shown is for 
the "across the grain" direction of the graphite which was the direc- 
tion of heat flow in the calorimeters. It can be seen in the figure 
that the experimental data used to determine the curve is scattered. 

The thermal conductivity of any graphite is difficult to control. 

The effects of this variation on the test data will be discussed in 
the "Accuracy of Results" section. The density and specific heat of 
graphites do not vary nearly as much and can be controlled readily. 

Rata Reduction 

All of the experimental data were recorded on magnetic tape so 
that the data could be reduced by computers. The thermocouple millivolt _ 
readings were converted into degree Fahrenheit according to National 
Bureau of Standards conversion tables. The temperature data were 


printed out in specified time increments over the firing time of the 
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rocket motors . The chamber pressure data were reduced in a similar 
manner . 

The heat flux into the calorimeter was determined by dividing 
the graphite core into 20 finite elements (fig. 13) , and specifying 
the temperature- time history of the third element. The temperature- 
time history was that of thermocouple 1, the thermocouple nearest the 
heated surface in each calorimeter. By using a finite element tech- 
nique, heat balance equations can be written for each of the elements 
over small time increments. For example, considering the heat flow 
to be one dimensional, the heat balance on element 1 is 


^in Sout + ^stored 


(16) 
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where is the temperature of element 1 calculated in the previous 

time interval. After the heat balance equations are written for each 
of the elements, these equations may be solved simultaneously for the . 
heat transfer coefficient. 

A computer program was used to solve for the heat transfer 
coefficient versus firing time for each of the calorimeters. The 
program utilizes the Gauss-Jordan method to solve the simultaneous 
linear heat balance equations . In addition to a temperature -time 
history of one of the elements, necessary input data include the 
physical dimensions of each element and the thermal properties of the 
material (K, p, and C ) which may vary with temperature. Output 
includes the heat transfer coefficient, the heating rate at the surface 



20 

of the calorimeter,, and the temperature gradient along the longi- 
tudinal axis of the calorimeter for each specified time increment. 

Radiative heating from the exhaust -gases was not accounted for in 
the analysis. This -was based on the fact, as discussed by Kuby (ll), 
that gases radiate energy in finite frequency bands as opposed to most 
solid bodies -which radiate in a energy continuum and therefore the 
total integrated value of emitted energy is much less for gaseous 
radiation. The propellant used in the tests was chosen because the 
amount of solid particles in its exhaust gases is small enough to be 
neglected. By not accounting for the radiation from the exhaust 
products, the data reduction is simplified, and the results are not 
effected substantially. 
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RESULTS 

Presentation and Discussion of Results 
Experimental heat transfer coefficients determined, for average 
chamber pressures of 220, 4l0, and 742 psia are presented. The experi- 
mental coefficients are compared 'with theoretical -coefficients deter- 
mined from the two methods of Bartz . In determining the theoretical 

\ 

coefficients with the simultaneous solution of the boundary layer 
momentum and energy equations, the skin friction coefficient was* 
evaluated at the arithmetic mean between the free stream and wall 
temperature as well as at the free stream temperature. Also, the 
theoretical coefficient was evaluated with an "interaction component" 
value of 0.1 as recommended by Bartz and also with a value of zero 
which essentially decouples the momentum and energy equations. 

Other data presented are the calculated temperature distribution 
derived from the temperature data from the thermocouple 1 of each 
calorimeter. This calculated temperature distribution is compared with 
the temperature readings of the thermocouple mounted along the longi- 
tudinal axis of the calorimeters to determine the accuracy of the 
experimental data. The deviations between thermocouple readings and 
the calculated temperatures are discussed in the "Accuracy of Results" 
section. The experimental heating rates are also presented. 

The thermocouples are numbered consecutively from one to five with 
thermocouple 1 located nearest the heated surface of each calorimeter. 


Average Chamber Pressure of 220 psia 

Experimental heat transfer data were obtained at a local area 
A 1 

to throat area ratio , -r— , of 1.785 in the convergent section of the 

•^t 
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‘ A 1 

nozzle, at the nozzle throat, -r— = 1, and In the divergent section at 
A 1 ** 

-t— = 2.569* Figures it through 22 show the temperature response of the 

. 

three thermocouples nearest the heated surface in each of the calori- 
meters, the experimental heating rates, and the experimental heat 
transfer coefficients. The oscillatory nature of the heating rate and 
heat transfer coefficient curves is due to the data reduction technique 
rather than physical conditions in the test. 

The heating rate curves show a decline in the magnitude of the 
heat flux with increasing time. This is to be expected since the 
driving potential, ^T^ - T^j , becomes smaller with increasing time. 
However, the heat transfer coefficients calculated from the measured 
temperatures in the convergent section and at the throat also decreased 
with time. In Figure 25 it can be seen that both of Bartz's methods 
overpredict the heat transfer coefficient throughout the nozzle. The 
heat transfer data from the latter part of the firing can be correlated 
by using a laminar flow equation used in (12) : 


N X If 0 ' 067 

St x iN Pr 5 Ke D V 


(18) 


This is Pohlhausen's equation for laminar flow where the characteristic 
length has been changed to the local diameter since the flow is 
internal. The heat transfer coefficient determined from the experi- 
mental data in the divergent region of the nozzle did not vary as much 
as the throat and convergent heat transfer coefficients. The predic- 
tions from the simultaneous solutions, where the skin friction 
coefficient was evaluated at the free stream temperature, and from the 
laminar flow equation bracketed the test data in the divergent region. 
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The leveling of the slope of the temperature curves of thermo- 
couple 1 in the convergent and throat calorimeters seemed unusual. In 
order to determine if the temperature data were valid, a temperature- 
time response of thermocouple 1 in each calorimeter was assumed as shown 
in Figures Zk and 25* These assumed temperature responses were used 
to compute the temperature response throughout the calorimeter with the 
use of the finite element computer program. Figures 2-4- and 25 shew that 
on comparing the results using the assumed response with -thermocouples 
2 and 5 in the convergent calorimeter and thermocouple 2 in the throat 
calorimeter, the assumed temperature response was incorrect. The 
calculated responses were higher than the experimental responses at 
the locations where temperature were measured. These comparisons 
lend credance to the experimental data which show that the heat transfer 
coefficient becomes smaller with time. 

It is known ( 4 ) that cooling the laminar boundary layer increases 
the stability, i.e., increases the critical Reynolds number for 
transition to turbulent flow. The boundary layer in this test was 
cooled by the nozzle wall throughout the test but to a lesser extent 
as time increased. Thus, it could be assumed that the conditions for 
laminar flow became less suitable at later times in the firing. The 
experimental data, however, indicate a tendency from turbulent or 
transition flow to laminar flow with' increasing, time. 

It is generally assumed ( 4 ) that boundary layer heat transfer coef- 
ficients are affected by wall temperature. This is accounted for by 
evaluating the skin friction coefficient and the gas transport properties 
at some intermediate temperature between the adiabatic and wall tempera- 
ture. Figure 25 shows the difference between evaluating the skin friction 
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coefficient at the free stream temperature and the arithmetic mean 
of the free stream and wall temperature for turbulent flow. Increasing 
wall temperature would tend to lower the heat transfer coefficient but 
not to the extent experienced in the test. 

A favorable pressure gradient tends to stabilize a laminar boundary 
layer and since the local pressure decreases along the length of the 
nozzle, a turbulent boundary layer could revert to laminar flow. In 
Figure 2 6 the experimental data from the 220 psia chamber pressure 
test is compared with experimental data presented in (12) . The figure 
is a plot of the parameter, Stanton number multiplied by the Prandtl 
number to the 0.6 power versus free stream Reynolds number based on the 
local diameter. It can be seen that the data reported on herein falls 
within the transition region as defined in (12) . These results suggest 
that it is possible for laminar or transition flow to occur in a rocket 
motor nozzle even though the exhaust gases originate from a combustive 
process within the motor chamber that is highly turbulent. 

Average Chamber Pressure of 4-10 psia 

The nozzle used in this test was the same as that used in the 
220 psia test with measurements made at. the same locations. The addi- 
tional pressure was generated by coupling two motor cases together 
with the head end case -containing only half as much propellant as the 
second case. This method gave a chamber pressure that varied more than 
in the 220 psia test. 

Figures 27 through 35 show the temperature response of the three 
thermocouples nearest the surface in each of the calorimeters, the 
heating rates, and the experimental heat transfer coefficients. 
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Figure 56 shows a comparison of the experimentally determined coef- 
ficients with the several methods of predicting the coefficients. It 
can be seen that, as in the 220 psia test, both of the Bartz methods 
overpredict the heat transfer coefficients, with the technique of 
evaluating the skin friction coefficient at the free stream temperature 
coming nearest to correlating the data. The heat transfer coefficient 
at the throat of the nozzle decreased sharply with time. This is also 
occurred to a lesser extent at the convergent calorimeter. The local 
Reynolds number at the convergent, throat, and divergent calorimeter 
was 0.725 X 10^, 1.02 x 10^, and 0.79^ X 10^, respectively, indicating 
turbulent flow when compared to the data in (12) . The drastic reduc- 
tion of the heat transfer coefficient at the throat is questionable 
in view of the fact that the Reynolds number at the throat indicates 
turbulent flow and that the magnitude of the coefficient drops slightly 
below that of the convergent calorimeter. The initial 1.5 to 2.0 seconds 
of data do indicate that the Bartz methods slightly overpredicted the 
results of the test at the throat. 

There was very little variation with time in the heat transfer 
coefficient on the divergent calorimeter. By evaluating the skin 
friction coefficient at the free stream temperature, the heat transfer 
coefficient was very nearly predicted. The convergent calorimeter 
data did show some decrease in the heat transfer coefficient. The 
heat transfer data at this point was slightly overpredicted by using 
the free stream skin friction coefficient in the simultaneous solution. 

As mentioned before, there was some variation in the chamber 
pressure. An average chamber pressure of 4l0 psia was used in the 
analytical solutions. The chamber pressure varied from 44o psia in the 
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initial stages of the firing to 5^0 psia just before burnout. Using 
the Nusselt number correlation equation, -where the heat transfer 
coefficient is proportional to the 0.8 power of the chamber pressure, 
to determine the effect on the heat transfer coefficient, it is found 
that the chamber pressure variation would theoretically cause the 
heat transfer coefficient to vary 6.5 percent above to 9 *8 percent 
below that for an average chamber pressure of 4l0 psia. The change 
in chamber pressure could partially explain the reduction in heat 
transfer coefficients at the convergent calorimeter and at the throat 
calorimeter. 


Average Chamber Pressure of lk-2 psia 

The nozzle contour for this test was slightly different for this 

test in order to increase the chamber pressure to 7^-2 psia. The 

contour still included at' ^5° half angle convergent region and a 15° 

half angle exit cone, but had a. smaller throat diameter . Heat transfer 

measurements were' made at an area ratio, — = 2. 067 in the convergent 
A, A/* ' • 

region, — - 1 , at the throat, and at -r— - 5 - 76 ^-, in the divergent 

^t ■ , \ 

region. ‘Figures 37 through ^5 show the temperature responses of the 


three thermocouples -nearest the heated surfaces in each of the calori- 
meters, the heat- fluxes calculated from the' temperature responses, and' 
the resulting heat transfer coefficients. .Figure k6 compares the experi- 
mental results with the various analytical prediction methods. This 
figure shows the data to be slightly overpredicted when the free stream 
skin friction coefficient is used in the simultaneous solution. The 
convergent region data showed very little variation of the heat transfer 
coefficient during the burning of the motor although the overprediction 
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■was the largest at this measurement location. The heat transfer coef- 
ficient at the throat was in the range of the predictions that used the 
free stream skin friction coefficient initially and decreased as the 
burning progressed. The divergent calorimeter data agreed very nearly 
with the predictions of the method stated above. 

Accuracy of Results 

The accuracy of the experimental results can be broken down into 
two parts: 

(1) Accuracy of the measured data. 

(2) Accuracy of the method of data reduction. 

Factors which affect the accuracy of the measured data are the ability 
to determine the true location of the thermocouples with respect to the 
nozzle internal wall and the error introduced by the recording equipment. 
Factors which affect the accuracy of the data reduction are how well 
the properties of the ZTA graphite can be determined, i.e., the density 
and the constant pressure specific heat and thermal conductivity as .a 
function temperature, and how accurate the method of data reduction is . 

The location of the thermocouples in the calorimeter could be 
determined to within,±.003 inch by using close tolerances in the fabri- 
cation >of the calorimeter and X-raying each calorimeter. Figure 47 is 
a typical X-ray photograph. The tolerances on the installation of 
the calorimeter in the nozzle were such that the true location of the 
thermocouples with respect to the nozzle internal wall could be deter- 
mined to within ±.008 inch. This resulted in a band of uncertainty of 
the temperature distribution along the length of the calorimeter as 
shown in Figure 48 . 
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The thermocouples used in t he tests could measure true temperature 
•within 1 percent. The noise in the recording channels of the measuring 
equipment could he filtered so that there was negligible error in the 
recording of the temperatures. 

As stated before, the experimental data on the thermal conductivity 
of ZTA graphite is scattered (Fig. 12). In attempting to fit the 
temperature data from the thermocouple 1 with the experimental 
temperature distribution, it was necessary to vary the thermal conduc- 
tivity of the graphite. The curves shown in the figures that present 
the temperature data versus time represent a "best fit." Since the 
thermal conductivity of the graphite varies with temperature, the 
thermal conductivity versus temperature curve had to be varied rather 
than a thermal conductivity value. The density and constant pressure 
specific heat versus temperature of ZTA graphite are accurately 
know (l8) and error in their values were not considered in the analysis. 

In order to determine the effect of the ±.008 inch error in the 
location of the thermocouple 1, an error analysis was made by 
varying its location in the data reduction. The calculated heat input 
to the nozzle wall was found to vary ±5 percent during the early times 
in the test with the variation decreasing to ±2 percent during the 
later times in the test. This variation was in the calculated heat 
input only. As previously stated, the experimental data was not 
perfectly fitted with the calculated data based on the temperature data 
of thermocouple 1. Figures 48 through 50 show the percent difference 
between the heat input based on the temperature data of thermocouples 
1, 2, and 5 and the calculated heat input for late times in the tests. 
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The percent difference was determined by comparing the slopes of 
curves at the heated surface of the calorimeter since the heat input is: 

( 19 ) 

surface 

The thermal conductivity was not varied since the difference- in surface 
temperature of the two curves was not large. 

In the figures that show the calculated heat input and the 
calculated heat transfer coefficients, there is oscillation in the 
curves caused by the data reducing process rather than any physical 
phenomena occurring in the tests. This was caused by the large change 
in temperature of the graphite in a short time, the temperature varia- 
tion of the specific heat and thermal conductivity of the graphite, and 
the fact that a finite difference technique was used. 

In summary, the heat transfer results from the tests were not 
exact results. However, the data did show trends which were presented 


• St 
" " K <5 


and discussed previously. 
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SUMMARY AND CONCLUSIONS 

The experimental heat transfer coefficients obtained by firing a 
solid propellant rocket motor at chamber pressures of 220, 410, and 
7^2 psia were generally below the predictions made by the Bartz methods 
of solving the boundary layer energy and momentum equations simultaneously 
and the Nusselt number. correlation equation. 

The test at 220 psia chamber pressure demonstrated that it is 
possible to obtain laminar or transitional boundary layer flow in a 
solid propellant motor even though a turbulent combustive process is 
occurring in the motor chamber. The best agreement was found in the 
divergent region of the nozzles at all three pressures where the experi- 
mental coefficients were only slightly below the predictions from the 
simultaneous solution and evaluating the skin friction coefficient at 
the local free stream temperature. 

Data from the convergent and throat calorimeters indicated the 
heat transfer coefficient decreases with increasing wall temperature. 

This variation was more pronounced than the variation in the predic- 
tions when the skin friction coefficient was evaluated at the arithmetic 
mean of the free stream and wall temperature (film temperature)". The 
predictions using the film temperature skin friction coefficient were 
always above the experimental heat transfer coefficients in magnitude. 
Predictions using the skin friction coefficient evaluated at the free 

stream temperature were in the same range of the experimental data, 

’ > 

but they do not account for wall temperature variation. 

These results tend to substantiate the data in the literature in 


the respect that predicting heat transfer rates in the inlet portion of 
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a rocket motor nozzle is difficult (4, 19, 20 ) and that transitional 
and laminar flow along the nozzle wall is possible (l, 9 , 12) . The 
data in (5> 12) were collected from nozzles similar to the nozzles 
used in the tests reported on here and for similar combustion chamber 
conditions (P^ on the order of 200 psia) . The data from (l) showed 
the transitional flow phenomena to occur at a lower combustion chamber 
pressure (75 psia) . 

These test data indicate that by using either the Bartz simultaneous 
solution or the Nusselt number correlation equation recommended by 
Bartz, calculations of the heating load along a nozzle wall would be 
conservative in that the analytical heat transfer coefficient is higher 
than the experimental coefficient. It is felt, however, that the degree 
of conservatism will hamper efforts to have an efficient nozzle design. 
Further analytical and experimental efforts are needed to refine and 
improve Bartz 1 s methods . 

It should be noted that these sets of data are not a complete test 
of the Bartz methods since the data was derived from the use of one 
type of nozzle. A complete conclusion may be drawn when test data are 
available from a variety of propellants and nozzles used on solid 


propellant motors. 



32 


LIST OF REFERENCES 


1. Back, L. H., Massier, P. F., and Gier, H. L. 1964. Convective 

Heat Transfer in a Convergent-Divergent Nozzle. Intern. J. 

Heat Mass Transfer 7(5) * 549-568. 

2. Bartz, D. R. 1955- An Approximate Solution of Compressible 

Turbulent Boundary-Layer Development and Convective Heat 
Transfer in Convergent-Divergent Nozzles. Trans ASME, 77(8-), 
1235-1245. 

3 . Bartz , D. R. 1957 • A Simple Equation for Rapid Estimation of 

Rocket Nozzle Convective Heat Transfer Coefficients. Jet 
Propulsion 27(l), 49-51. 

4. Bartz, D. R. 1965 . Turbulent Boundary-Layer Heat Transfer From 

Rapidly Accelerating Flow of Rocket Combustion Gases and of 
Heated Air. Advances in Heat Transfer, Vol. 2. Edited by 
James P. Hartnett and Thomas F. Irvine, Jr., Academic Press, 

New York. 

5. Brinsmade, A. F., and Desmon, L. G. 1964. Hypothesis for Correlat- 

ing Rocket Nozzle Throat Convective Heat Transfer. Chemical 
Engineering Progress Symposium Series, No. 59* Vol. 6 l, 88 - 98 . 

6. Brokaw, R. S. 1961. Alignment Charts for Transport Properties 

Viscosity, Thermal Conductivity, and Diffusion Coefficients 
for Nonpolar Gases and Gas Mixtures at Low Density. NASA 
Technical Report R-8l. 

7 . Coles, D. E. 1961 . The Turbulent Boundary Layer in a Compressible 

Fluid. Rand Corporation, Report Number P-2417, Santa Monica, 
California . 

8. Elliot, David G., Bartz, Donald R., and Silver, Sidney. 1963 . 

Calculations of Turbulent Boundary-Layer Growth and Heat 
Transfer in Axisymmetric Nozzles. Jet Propulsion Laboratory. 
Technical Report No. 32-387* Pasadena, California. 

9 . Fortini, A., and Ehlers, R. C. 1963 . Comparison of Experimental 

to Predicted Heat Transfer in a Bell-Shaped Nozzle With 
Upstream Flow Disturbances. NASA TN D-1473* 

10. Kolozsi, J. J. 1958. An Investigation of Heat Transfer Through the 

Turbulent Boundary-Layer in an Axially Symmetric, Convergent- 
Divergent Nozzle. Masters Thesis, Dept, of Aeronautical and 
Astronautical Engineering, Ohio State University, Columbus, 
Ohio. 

11. Kuby, W. C. 1964. The Internal Environment of a Solid Propellant 

Rocket Nozzle. AIAA Preprint No. 64-158. 



33 


12 . Lee, Roland E. 19^3 • Heat Transfer to the Throat Region of a 

Solid Propellant Rocket Nozzle. United States Naval Ordnance 
Laboratory, TR 62-72, White Oak, Maryland. 

13* McAdams, W. H. 195^ • Heat Transmission, 3rd ed. McGraw-Hill 
Book Company, New York. 

l4. Shapiro, A. H. 195^ « The Dynamic and Thermodynamics of Compressible 
Fluid Plow, Vol. 2. The Ronald Press Company, New York. 

15- Sibulkin, M. 1956. Heat Transfer to an Incompressible Turbulent 
Boundary Layer and Estimation of Heat Transfer Coefficients 
at Supersonic Nozzle Throats. J. Aeron. Sci. 23(2), 162-172. 

16. Svehla, R. A. 1962. Estimated Viscosities and Thermal Conduc- 

tivities of Gases at High Temperatures. NASA Technical 
Report R-132. 

17. Union Carbide Corporation. 196^. Industrial Graphite Engineering 

Handbook. Union Carbide Corporation, New York. 

18 . Wagner , P . , and Dauelsberg , L . B . 1966 . Thermal Conductivity of 

ZTA Graphite. Los Alamos Scientific Laboratory. Report 
Number DC-8l63> Los Alamos, New Mexico. 

19. Welsh, W. E., Jr., and Witte, A. B. 1962. A Comparison of Analyti- 

cal and Experimental Local Heat Hexes in Liquid-Propellant 
Rocket Thrust Chambers. J. Heat Transfer c84-'(l), 19-28. 

20. Witte, A. B., and Harper, E. Y. 1963* Experimental Investigation 

of Heat Transfer Rates in Rocket Thrust Chambers. AIAA J. 

1 ( 2 ), kk3-k51. 

21. Zeleznik, Frank J., and Gordon, Sanford. 1962. A General IBM 70^ 

or 7090 Computer Program for Computation of Chemical Equili- 
brium Composition, Rocket Performance, and Chapman- Jouguet 
Detonations, NASA TN D-l4-5^ . 



33 ^ 


APPENDIX 




to Achieve Different Chamber Pressures 


Combustion , Pounds of 
Chamber Volume, in ? Propellant 

322 6.75 
6bk 10.13 
322 6.75 


- ( - 



35 


Table 2. Exhaust Gas Properties 


Average chamber 
pressure, psia 

Chamber gas 
temperature, °R 

Specific heat 
ratio 

Specific heat 
constant. Pressure, 
Btu 

lb °R 

. . lb 
Viscosity, 

Thermal conductivity 
Btu ft 
ft 2 sec °R 


220 

4io 

742 

4845 

488o 

4900 

1.23 

I.23 

I.23 

0.¥i3 

.0.443 

0.443 

48.6 x 10“ 6 

48.6 x 10“ 6 

48.6 x lo -6 

35-8 x 10" 6 

35-8 x 10" 6 

35.8 x 10“ 6 


Exhaust gas 
constituents,, mole 
fractions 


CO 0.193 * 

CL 0.005 ’ 

C0 2 4). 082 * 

H O'. 005 

Hq 0 .119 * 

HCL O.I 67 * 

HoO O .356 * 

1*2 0.089 * 

OH 0.004 


Note: 


Mole fractions of constituents did not vary with chamber pressure 




(a) Control volume of real flow (t>) Control volume of potential flow 


Figure 1. Control volumes for real and potential flow in nozzle 



Dimensions In Inches 


Floy Conditions 

P o = 300 psia C p = 0.567 Btu/lb °R 7=1-2 

3? 0 = t-,500° R K = O.83 \i ~ 1 JP ’°5 

p 0 = 1.3 X 10 - 8 lb sec/ft 2 

Figure 2 , Nozzle for evaluation of constant in heat-transfer coeffi 
cient equation (from ref. !{•) 








Chamber pressure, psia 







Figure 7 • Calorimeter locations within nozzle 




















Note: All thermocouples -were located on graphite core centerline 
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Five thermocouples, 
equally spaced 0.28 apart 


Five thermocouples 
equally spaced 0.233 apart 


Dimensions for 
convergent and throat 
calorimeters 





Dimensions for 

divergent 

calorimeters 

Silica phenolic 
sleeve -with 
0.10 wall 

ZTA graphite 
core 
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All dimensions in inches 


Figure 8. Calorimeter dimensions 






Figure 9. Experimental solid-propellant motor mounted in thrust stand 
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Figure 10. Instrumented nozzle on experimental solid-propellant motor 
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Convergent and throat 
calorimeters 


Element 

Thickness 

' 1 

.0210 

2 

.0210 

3 

.0360 

4 

.0475 

5 

.0623 

6 

.0625 

7 

.0625 

8 

.0625 

9 

,0625 

10 

.0625 

11 

.0625 

12 

.0310 

13 

.0625 

l4 

.0940 

15 

.0625 

1 6 

.0625 

17 

.0623 

18 

.1040 

19 

.1040 

20 

•lo4o 


Divergent 

calorimeters 


Heat in 



Side and end 
perfectly insulated 


Lement 

Thickness, in 

1 

.0210 

2 

.0210 

3 

.0360 

4 

.0475 

5 

.0625 

6 

.0625 

7 

.0860 

8 

.0388 

9 

.0623 

10 

.0625 

11 

.0625 

12 

.0625 

13 

.0625 

l4 

.0625 

15 

.0623 

16 

.0623 

17 

.0623 

id 

.0314 

19 

.0314 

20 

.0314 


Figure 13, ZTA graphite core finite elements 
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Figure 14. Comparison of experimental and calculated temperature response of thermocouples 1, 
2, and 3- Convergent calorimeter - P c = 220 psia 
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Figure 15. Heating rate vs burn time. Convergent calorimeter - P Q - 220 psia 
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Figure 17 . Comparison of experimental and calculated temperature response of thermocouples 1 
and 2. Throat calorimeter - P c = 220 psia 
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Figure 20. Comparison of experimental and calculated temperature response of thermocouples 1, 
2, and 3. Divergent calorimeter - P c = 220 psia 
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Figure 21. Heating rate vs burn time. Divergent calorimeter - P c = 220 psia 
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Figure 22. Convective heat -transfer coefficient vs burn time. Divergent calorimeter 
P c = 220 psia 
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Figure 2b. Extrapolated thermocouple number 1 response. Convergent calorimeter - P c = 220 psia 
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Figure 25- Extrapolated thermocouple number 1 response. Throat calorimeter - P = 220 psia 
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Figure 27. Comparison of experimental and calculated temperature response of thermocouples 1, o^ 

2, and Convergent calorimeter - P c = J+10 psia ^ 
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Figure 30. Comparison of experimental and calculated temperature response of thermocouples 1 , 
2 , and 3- Throat calorimeter - P c = 410 psia 
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Figure 31. Heating rate vs burn time. Throat calorimeter - P c = 410 psia 
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Figure 32. Convective heat -transfer coefficient vs burn time. Throat calorimeter - _ct\ 

P c = klO psia 
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Figure 35- Convective heat -transfer coefficient vs burn time. Divergent calorimeter 
P c = 410 psia 
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Figure 37 * Comparison of experimental and calculated temperature response of thermocouples 1 
2 , and 3* Convergent calorimeter - P„ = 74-2 psia 
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Figure 38. Heating rate vs "burn time. Convergent calorimeter - P c = "jkZ psia 
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Figure 39* Convective heat-transfer coefficient vs burn time. Convergent calorimeter - 

P c = 742 psia ^ 




Figure 40. Comparison of experimental and calculated temperature response of thermocouples 1, 
2 , and 3* Throat calorimeter - P c = 7^-2 psia 
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Figure 42. Convective heat-transfer coefficient vs burn time. Throat calorimeter 
P c = 742 psia 
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Figure 43- Comparison of experimental and calculated temperature response of thermocouples 1, 
2 , and 3. Divergent calorimeter - P c = 742 psia 
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Figure 45- Convective heat-transfer coefficient vs hum time. Divergent calorimeter - 
P c = "jk-2 psia 
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Figure k 6 . Comparison of ejqoerimental and analytical heat-transfer coefficients - p„ = 742 psia 
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Figure 47. Typical X-ray of calorimeters 
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Figure 50. Comparison of analytical and experimental temperature distributions - P c = Jb-2 psxa 


